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' We have investigated the effect of weak disorder on the ground state of a two dimensional electron 

gas in the quantum Hall regime at filling factors slightly deviating from unity. The skyrmions and 
antiskyrmions are found to be present even at filling factor v = 1. They may be strongly squeezed 
by the disorder. We have studied temperature effect on skyrmion density and size. 
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£C) \ Recently it has been proved that the lowest energy charged excitations of a two-dimensional electron gas (2DEG) 
in quantum Hall regime at filling factor v — 1 are skyrmions 0] . These excitations involve many electrons and can be 
|_i ' adequately described in the framework of a semiclassical cr-model as the topologically non-trivial spin textures of an 
isotropic ferromagnet ||] . Owing to the incompressibility of the Quantum Hall ground state the topological charge of 
r*| | a skyrmion is equal to its electric charge || . 

There is a convincing experimental evidence for these collective excitations. However, the experiments show 
only qualitative agreement with the theory ||||. It has been suggested |||| that the agreement will be better if the 
effects of disorder are taken into account. Despite the importance of the topic, very few theoretical studies [|]|| deal 
with skyrmions in disordered heterostuctures. The weak disorder limit has received no attention. 

In this letter we show that even a weak disorder gives rise to (anti)skyrmions even at filling factor v = 1. In the 
limit of a small Zeeman energy the Skyrmion radius is decreased by the disorder potential. We present the results for 
I . the ( ant i) skyrmion density at filling factors deviating from unity. We have found an interesting temperature effect on 
skyrmion density and radius. We attempted to compare our results with experimental data of M. 



C 

O 

o 



In the limit of small g, which stands for the ratio between the Zeeman and exchange energy, skyrmions can be 
described with an effective energy functional derived in . The energy of a single skyrmion as a function of its radius 
a reads 
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where Eq is the skyrmion gap energy, Is is the magnetic length, g is the electron g factor, r s ~ Ib!) 1 ^ 2 ■ Although the 
gap energy dominates the total energy, the radius is set by the two smaller terms depending on a. The equilibrium 
Skyrmion radius ao can be found by minimizing the total energy with respect to a. In this article we neglect the 
Qv^ ' slowly changing logarithmic factor in the Zeeman energy and incorporate it to g. 

The disorder is brought into a 2DEG mainly by the donors situated in a layer several hundreds of angstroms away 
from the electron layer. We incorporate the effect of disorder in a fourth term in the energy functional. For a Skyrmion 
at position ro in a random donor potential V(r) this energy reads 



Sdis(ro) = f d 2 rp{v - r , a)V(r) (2) 

where p(r, a) = (a/r 2 + a 2 ) 2 /ir is the charge density of a Skyrmion with radius a. The potential V(r) is not screened 
by the electrons of the 2DEG, since at integer filling factors the 2DEG is an insulator. We assume the Gaussian 
distribution of V(r) with a correlator 
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V q f) = f (3) 

which is calculated by averaging over donor positions. Weak disorder requires Uq <C Eq. If we assume that donor 
positions do not correlate, we come up with an unrealistically big estimate of Uq ~ e 2 fed, d being the distance 
between the 2DEG and the donor plane. In realistic heterostructures the observed disorder is known to be much 
weaker, possibly due to correlations between donor positions. Therefore, we assume that Uq <C e 2 /d, Eq. Uq can be 
extracted from the experimentally measured mobility at zero magnetic field. 

At a very qualitative level, our results are straightforward. At zero temperature it is energetically favorable for a 
Skyrmion to exist in the potential wells which are sufficiently deep (Fig. 1), so that -Edis < Eq. Since the disorder is 
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weak, such wells are infrequent, so that the Skyrmion density will be much smaller than the electron density. These 
deep wells are to be expected very steep, that compresses the skyrmion to make it fit into the well, (see Fig. 1) This 
agrees with the preliminary results of computer simulations reported in || . 

At temperatures of the order of Uq the skyrmions may appear in the wells which are less deep (left side of Fig.l), 
that increases their density. Since the potential minima become less steep, the skyrmion radius becomes bigger. The 
same consideration is valid also antiskyrmions which are situated in the maxima of the random potential. 

Let us make these simple ideas quantitative. To calculate the skyrmion density at filling factor v = 1 and zero 
temperature, we consider the following integral 

n sk {E, a)= f d 2 r £ <5(r - r m )6(a - a m )S(E dis {r, a) - E)/ f d 2 r (4) 

where n sk (E,a) is the density of skyrmions having potential -Edis = E with radius a per interval of energy and 
radius (that we can call skyrmion density of states) and {r m ,et m } is a set of Skyrmion positions and radii for a 
certain V r (r)-configuration. These {r m ,a m } are not arbitrary, but determined from the fact that the skyrmions are 
in equilibrium. Since we assume that Skyrmion density is small, we can consider isolated Skyrmions. It means that 
we disregard interactions between Skyrmions when they are in different wells, but due to strong repulsive interaction 
within the well there is only one Skyrmion in each well. In this case a Skyrmion only interacts with the disorder 
potential. This results in three conditions, two for Skyrmion position: d x Edi s = dyE^is = and one for Skyrmion 
radius: 9 a ^dis = —d a E sk . The actual density of skyrmions is to be found by integrating n sk (E, a) over all radii and 
E < E sk rts Eq. 

To proceed we rewrite this integral as follows. To incorporate the equilibrium conditions we mentioned above, 
we transform (5-functions in terms of the three equilibrium conditions. The Jacobian of this transformation enters 
the integral: \d a p(E s k + Edi S )\, with a,/3 — {x,y,a}. Them we average Eq. (Q) at r = over all possible V q - 
configurations. We have to take into account that skyrmions are situated in minima, so that that all eigenvalues of 
the Jacobian are positive. We express E^ in terms of Fourier components: -Edis(O) = J2 ^dis(q), with Ed is (q) = V q p q 
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and p q — aqKi (aq) By adding extra variables we can rewrite all 5- and ^-functions in the form of exponents. 

Then integration over the random potential appears to be Gaussian and can be performed. In the course of averaging 
the terms like ^2 p 2 < V q 2 > will appear. Formally they diverge at q — > 0. To deal with this divergence we have 
i 

to recall that the interaction between Skyrmions is a long range one. So that it becomes effectively important at 
large distances and Skyrmions screen out components of disorder potential having q < ^/n sk , inverse distance between 
Skyrmions. So we cut off integration at q ~ ^Jn sk . 

After integration over extra variables the expression for density of states reads 
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where 



a = 2tt ^ ~ - ln(a^/nZ) (6) 
q 

Here we take into account that the sum in Eq. || diverges and must be cut off at q ~ ^/n sk . We assume that 
ay/n sk <C 1. So that a depends on Skyrmion concentration and eqs.(||),(^) form together a self consistency problem. 

Now we evaluate the radius a which minimizes the exponent of eq. (q). We substitute the Skyrmion energy of 
eq. (|l|) in the exponent of eq. (||) . It appears that there are two limits corresponding to the strength of the Zeeman 
energy 
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where ao is the skyrmion radius in the absence of a random potential field. 



2 



Now we can proceed with calculation of the total density. We expand the exponent around a — a max and E = —E s ^. 
We restrict ourselves to the more interesting limit of small Zeeman energy: <7 3 |j^ *C 1. In this limit we encounter 
again two different limits depending on the strength of the Zeeman energy. The point is that the prefactor, which 
is in fact the determinant \d a p(E s ^ + E^ s )\, vanishes at a — a opt in the first order approximation. In both cases 
the exponent is the same, n s k oc exp(2-7ri?Q/a). This is sufficient to solve the self consistency problem. The leading 
approximation to a is ctQ = ^/ttj^, the next order approximation is necessary to evaluate the prefactor. Finally we 
obtain for the density 
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where n e is the electron density in the 2DEG, 71 ~ 9 x 10 -4 , 72 — 6 x 10 3 are numerical factors. Here we have taken 
into account that the number of skyrmions at v = 1 is equal to the number of antiskyrmions, and those two equally 
contribute to n s k- 

Therefore we find three regimes in dependence of the relative strength of the Zeeman energy with respect to disorder. 

The most probable Skyrmion radius is equal to its disorder less value at g^ jjft 3> 1 and to a opt — ^j-\J ^j~^b <C clq in 
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the opposite limit. At g§ jjft — f^H 3 tne prefactor changes. 

Let us evaluate skyrmion densities at filling factors slightly deviating from unity. In the absence of disorder, 
skyrmion and antiskyrmion densities are proportional to the filling factor deviation (Fig. 2), n s k,ask = ±Sf6(±df). 
where Sf = n e {y — 1). This is no longer correct in the presence of disorder potential. We note that at unity filling 
factor the chemical potential lies precisely in the middle of the gap between skyrmion and antiskyrmion states. If the 
filling factor deviates from unity, we have to recalculate the densities for a shifted chemical potential [i. Fortunately, 
this is simple. The density of states exhibits exponential dependence on energy, n s k(E) oc exp(— E/d), = 2oUq/Eq. 
This is why the skyrmion density is changed by a factor exp(— n/fl) whereas the antiskyrmion one by exp(/i/f2). An 
attention should be paid to the fact that the change of the chemical potential also affects a. We note that the total 
charge per unit area must correspond to the filling factor, that is — n as k = Sf. This allows as to express (i in 
terms of Sf and obtain for the densities 
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These densities are plotted in Fig. 2, together with the ratio of the total number of spin flips n s u + n as k and the 
number of spin flips at Sf = 0. The latter quantity can be observed experimentally since it determines the rounding 
of the spin polarization peak. || 

Let us consider the effect of a finite temperature. Still we have to assume that the Skyrmion density is much smaller 
than the electron density, that is, T <C Eq. Also the temperature should be low enough for Coulomb interaction to 
make it impossible for two Skyrmions to occupy the same potential well. Under these circumstance the skyrmions 
effectively behave as fermions, so we can make use of Fermi statistics. To calculate the total density at \i — we 
integrate the density of states of Eq. @ multiplied by the Fermi distribution function, risk(E, a) x f(E — i?o)(see Fig. 
3 a,b). If we expand the exponent in n s -^{E") near E = E^ yielding a term linear in E — E^ we see that the integral 
converges only at T < T C ,T C = Uo/i^/n. 

Of course, it does not mean that the density actually diverges at T = T c . Rather, it indicates an exponential 
temperature dependence of the density at T > T c . To reveal this one, we go back to the expression (||) and consider 
the most important term in the exponent which is quadratic in energy. This term competes with the exponent of the 
Fermi distrubution tail. The minimum is achieved at E = E m i n = aT^T . The solving the self consistency problem 
gives a(T) = E /(2T + 2T C 2 /T). At the critical temperature E min = E , as it should be. The density at T > T c thus 
reads 

n(T) cx eM-Eo T2 ^ T2 ) (10) 

and exhibits a non-Arrenius behavior. The logarithm of the density is plotted in Fig. 3 versus inverse temperature. 
The optimal size of a ( ant i) skyrmion can be extracted from the exponent (|^). At T > T c and it grows with temperature: 
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a(T) — a(0)T/T c , the latter is valid till a(T) <C do- This is because the potential minima become less steep at higher 
energies. 

With our method, we can also obtain the finite temperature results for v ^ 1. 

Our results can be checked with spin polarization measurements. Indeed, the finite density of (anti)skyrmions at 
v = 1 would manifest itself as a rounding of the spin polarization peak. The reduction of their size a due to disorder 
and its restoration at T > T c can be detected as a change of skyrmion spin oc a 2 . However, this check requires a very 
accurate measurement of the spin polarization peak in the close vicinity of v = 1 at different temperatures, that has 
not been performed yet. 

As to transport measurements, we attempted to compare our results with experimental data of the Ref. ||. We 
assume the longitudinal resistivity to be proportional to the skyrmion density: p xx ~ n(T) and the temperature 
dependence of n(T) to be dominated by the exponent ([h]). In Fig. 3 the experimental data are fit to the theoretical 
curve, with fit parameters Uq ~ 15 K and Eq ~ 12 K. Although Uq ~ Eq, we expect our theory to be qualitatively 
true since the numerical factors provide big exponents. Indeed, in this case T c ~ 2K <C Eq. From the given mobility 
we extract Uq ~ 10K. 

However, most of the transport measurements exhibit no saturation of p xx down to 1 K. Although this can be 
explained by better quality of the heterostructures (Uq < 3K), we hesitate to make a point out of our fit. The point is 
that all transport measurements, whatever interpreted, give an estimation of Eq which is an order of magnitude smaller 
than the theoretical value. Possibly this indicates that the assumption p xx ~ n(T) is no good. Our results explain 
how this may happen. Indeed, even at zero temperature and v = 1 we have a finite concentration of (anti)skyrmions. 
What may be thermally activated is their mobility rather than concentration. It is no surprise since they have to 
overcome high potential barriers when moving between potential minima. 

In conclusion, we have developed the theory of skyrmions in weakly disordered heterostructures. The skyrmions 
appeared to be present at filling factor v = 1. The disorder and temperature strongly affect their density and size. 
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FIG. 1. Skyrmions in disorder potential. Their size must be small to fit into a deep well. 



FIG. 2. Skyrmion densities versus filling factor. 



FIG. 3. Temperature dependence of the skyrmion density. Solid line presents our theoretical results, the squares correspond 
to the fitted data of Ref. 6. 
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